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HEAT AND MOMENTUM TRANSFER BETWEEN A SPHERICAL 
PARTICLE AND AIR STREAMS 


By Y. S. Tang, J. M. Duncan, and H. E. Schweyer 
SUMMARY 


Heat-transfer coefficients for a spherical particle heated by an 
induction coil in a moving air stream were experimentally determined 
for the Reynolds number range from 50 to 1000 using spheres of 1/8- to 
5/8-inch diameter and air velocities from 1 to 13 feet per second. A 
correlation of the heat-transfer factor or Stanton number with the 
Reynolds number was obtained and expressed by an empirical equation. 
This correlation is in agreement with the values calculated from theory 
for the lower range of Reynolds numbers studied. 


The skin-friction factor representing the momentum transfer calcu- 
lated from the boundary-layer theory shows good agreement with the experi- 
mental heat-transfer factor except in the lower range of Reynolds numbers 
studied. The relationship St = Cg/2 where St is the Stanton number 
and Cr is the skin-friction factor is suggested for the case of an air 
stream flowing around a sphere. 


An empirical equation relating the heat-transfer factor to the total- 
drag coefficient is also suggested. 


INTRODUCTION 


For several years, there has been considerable theoretical and 
experimental interest in heat transfer and momentum transfer (fluid fric- 
tion) for bodies submerged in a flowing fluid. These transfers occur 
frequently in engineering operations. They are becoming of increased 
importance in catalytic operations, flow in packed beds, calcining, gas 
absorption, combustion chambers, and other solid-gas and liquid-gas 
reactions. 


Johnstone, Pigford, and Chapin have made an analytical study of the 
heat transfer between a small spherical particle and ambient fluid stream 
(reference 1). Drake, Sauer, and Schaaf recently have made the same theo- 
retical analysis following Johnstone's assumptions but using a different 
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method for solving the differential equation (reference 2). Their 
solution, which differs from Johnstone's for Reynolds numbers based on 
the particle diameter Re below 1000, was taken as the theoretical 
basis for the present investigation. An extensive survey of experi- 
mental data in the literature has been made by Williams (reference 3) 

and an empirical curve correlating Nu with Re was recommended for 

a spherical particle in an air stream. However, discrepancies between 
the curve and the data are quite large in the region where the Reynolds 
number Re is less than 1000. Furthermore, those data used were either 
reported by the original observers as unreliable in this range or the 
Observations were made on a water system. Because of the difference in 
Prandtl number of the two fluids the results would not be expected to 
fall on the same curve (see reference 4). This region where the Reynolds 
number Re is less than 1000 was, therefore, selected for further study. 


The similarities between the heat- and momentum-transfer processes 
at the interface of a solid and fluid have been known since Reynolds 
first suggested the analogy theory (reference 5). This theory has been 
modified by various investigators, incliding Colburn (references 6 
and 7), Prandtl (reference 8), Taylor (reference 9), Von Kerman (refer- 
ences 10 to 12), and, more recently, Martinelli (references 13 and 14) 
‘and Boelter, Martinelli, and Jonassen (reference 15). It has been 
reported as a satisfactory principle in dealing with fluids flowing 
parallel to surfaces such as flat plates, conduits, or other confining 
surfaces (references 6, 14, and 16). Few comparisons had been made in 
the case of blunt objects submerged in fluids until the recent work on 
flow around cylinders (references 17 to 19). The difficulty in a study 
involving a blunt object lies in the separation of the effect of the skin 
friction from the total drag force exerted on the body. Sherwood, in 
his paper reviewing the relationship of these transfers in turbulent 
flow, compared the calculated friction data with heat-transfer data for 
flow around single cylinders and showed good agreement at certain Reynolds 
numbers (reference 20). However, since such a comparison for spheres had 
not been made previous to this investigation as far as is known, the skin 
friction on the surface of a sphere is calculated, and the skin-friction 
factor Cr is compared with the heat-transfer factor St. 


Since the total-drag coefficient is measurable and is available in 
the literature (references 21 to 23), a correlation between the total- 
drag coefficient Cp and the heat-transfer factor St would be of 
interest. This relationsip is developed for the range of Reynolds 
numbers Re studied where the turbulence at the wake of the sphere is 
a minor factor in the total drag. 


This investigation was conducted by the Department of Chemical 
Engineering of the University of Florida under the sponsorship and with 
the financial assistance of the National Advisory Committee for 
Aeronautics. 


NACA TN 2867 


SYMBOLS 


The symbols used for general correlations are dimensionless and 


may be in any consistent system of units. The experimental data are 
all expressed in English units and so are the computations herewith. 
All the symbols used are listed as follows: 


A 


Asep 


Ag,Ao 


B,C 


B(g) ,Βα(α) 


Kg 


kp 


area 

surface area up to separation point 

cross-sectional area and surface area, respectively 
constants in equation (Ak) 


functions of g acording to equations (BT) and (B8), 
respectively 


constant in equation (A2) 


total-drag coefficient [159081 drag force 
(ovo*] 2)As 


skin-friction factor Total frictional force 
GUI 


momentum-transfer factor 

Specific heat &t constant pressure 

diameter of spherical particle 

variable defined as πεαθ/ Εξ 

heat-transfer coefficient 

heat-transfer coefficient due to free convection 
plate current of oscillator 

thermal conductivity 

thermal conductivity of steel 


thermal conductivity of fluid 


Nu 


Po 


Tair 
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length of wire 

finite section of length Ρ 

length of volume element in y-direction 

constant 

constant in equation (2) 

Nusselt number  (hD/k) 

pressure 

static pressure in undisturbed stream 

perimeter 

Prandtl number (σρη| κ) 

rate of total heat transferred, which is equal A rate of 
heat generation in particle in equilibrium state 

total heat generated per unit time in sphere with radius Ro 

rate of heat transferred by conduction 

rate of heat transferred per unit area 


rate of heat transferred by radistion per unit area 


radius 


spheric coordinates 


perpendicular distance from axis of revolution to surface 
of body of revolution 


outside radius of spherical particle 
Reynolds number based on diameter of particle  (pDUg/L) 
heat-transfer factor or Stanton number (h/CpPUo) 


absolute temperature; in appendix A, temperature of wire 
referred to surroundings 


temperature of air stream 
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T1;To 
If,tg 
Tg 


η Ὅρ 


temperatures at stations 1 and 2, the hot and cold junctions, 
respectively 


film temperature and fluid temperature, respectively 
temperature at surface of coated thermocouple wires 


initial temperature and average temperature of particle, 
respectively 


temperature at center and at surface of particle, 
respectively 


temperature inside and outside of particle, respectively 
temperature of fluid outside the thermal boundary layer 
temperature difference between particle and main air stream 
velocity at outside of boundary layer 


average velocity of undisturbed stream 


maximum velocity at center of column 


average velocity over entire cross section of flow 


velocity of fluid in direction of spheric coordinates r, 
θ, and Ø, respectively 


component of local velocity at any point in x-direction 
component of local velocity at any point in y-direction 
plate voltage of oscillator | 
distance parallel to surface, orthogonal coordinate 
distance perpendicular to surface, orthogonal coordinate 
thermal diffusivity (κ/Ορο) 

hydrodynamic boundary-layer thickness 

thermal boundary-layer thickness 


time 
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H absolute viscosity 

ν kinematic viscosity | s 
p density 

T skin friction or shearing stress 


EEAT-TRANSFER APPARATUS 


Measurements of heat transfer were made on carbon-steel spherical 
particles suspended in an air stream and heated by an induction coil. 
The arrangement of the apparatus for this study is shown by figures 1 
and 2. Three particle sizes (1/8-, 1/h-, and 5/8-in. diameters) were 
used. A single particle was suspended by means of very fine thermo- 
couple wires in a vertical insulated glass column through which a stream 
of air was blown. A high-frequency (300-kc) induction coil which was 
wound around a section of the column was used to generate the heat in 
the particle. Two sizes of glass columns (2- and 3-in. inside diameters) 
were used so that the value of the Reynolds number Re could be changed 
either by changing the size of the particle or by changing the cross 
section of flow while the flow rate remained within the range of the 
rotameter. A change of the flow pattern of the main stream at a single 
value of Re could also be obtained by changing the column size. 


The high-frequency alternating current was generated by an oscil- 
lator (Westinghouse R-F Generator). Values of the plate voltage and 
current of the oscillator were used as a measurement of the energy input 
to the heating unit. The use of an induction coil as a heat source for 
the study of heat transfer was illustrated by Kramers (reference 24). 
This method was found advantageous because it allowed the use of an 
extremely small particle and avoided the difficulties that may be 
encountered with direct electric heating. 


The flow of air was obtained from a centrifugal blower discharging 
into a surge tank where a constant temperature (15° F) of the air was 
maintained by a thermostat and entered into the glass column through a 
Wire screen. The average velocity of the air stream in the column was 
measured by a calibrated rotameter and the flow rate was regulated, for 
each test, by a gate valve downstream from the rotameter. 


The actual velocity at the center of the glass column was measured 
by means of a pitot tube and a diaphragm-type micromanometer (reference 25) 
which proved to be more satisfactory than a tilting-type inclined dif- 
ferential manometer (reference 26) which was initially used. 
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The temperature of the particle relative to the air temperature was 
measured by a thermocouple which was also used to support the particle, 
with its hot junction located inside the particle and its cold junction, 
in the air stream. Number 40 gage Cupron (copper, 55 percent; nickel, 

45 percent) and Cromin D (nickel, 30 percent; chromium, 5 percent; and 
iron, 65 percent) wires were used in order to minimize the wire effect 

at the wake of the particle and the heat loss through conduction. Both 
thermocouple wires were insulated by coating with shellac. The potential 
difference of the thermocouple was determined by a Rubicon Precision 
Potentiometer. The junctions of the thermocouple were formed by elec- 
trical welding using a 6-volt direct-current source as recommended by 
Carbon, Kutsch, and Hawkins (reference 27). The hot junction was attached 
inside the particle near the center so as to have the least possible dis- 
turbance of the spheric surface. Several methods of attaching the junc- 
tion were tried. The method selected as most convenient was to drill 
and tap a hole to the center of the particle. The thermocouple wires 
were then passed through a steel screw, knotted on its end, and pushed 
to the center of the particle by the screw which was used to plug the 
hole in the particle. A sectional view of the particle with the thermo- 
couple junction attached is shown in figure 3. To insure a good contact 
between the body of the particle and the junction, a small drop of mercury 
was placed in the hole before the plug was inserted. This method has 

the advantage of making possible the use of the same thermocouple for 

the different particles except the smallest particle size used for which 
direct soldering of the junction to the inside of the particle was 
employed because of mechanical difficulties in applying the above method. 


Two copper-constantan thermocouples were uséd to obtain the tempera- 
ture of the air stream for making temperature corrections and evaluating 
properties of the fluid. One was placed in the test section of the 
column and the other at the outlet of the rotameter. 


TEST PROCEDURE 


The spherical particle was first placed in position inside the glass 
column and air was blown through the column with the flow rate controlled 
to a desired value which was measured by a carefully calibrated rotameter. 
After the flow reached a stable condition, current was applied to the 
heating unit. A desired constant energy input to the heating unit was 
obtained by adjusting the input resistance. The temperature difference 
was read at different intervals of time. The temperature of the 'par- 
ticle increased until the rate of heat loss from it was equal to the 
rate of heat generation in the particle. About 30. to ΠΟ minutes were 
usually required befcre this steady state of heat flow was8 reached. The 
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rate of heat transfer from the particle to the air stream was then 
determined from the rate of heat generation in the particle which is 
related to the energy input to the heating unit by an experimentally 
determined calibration (see section entitled "Preliminary Measurements" ). 
Because of the small temperature difference between the particle and 
surroundings it was considered permissible to assume the radiation loss 
from the particle to be negligible in comparison with the convective 
transfer. However, in order to evaluate any possible errors, calcula- 
tions were made for the radiation loss during cooling of the particle 
While it was enclosed in an evacuated glass bulb. These calculations 
and also a consideration of possible conduction loss through the sup- 
porting wires are treated in detail in the section entitled "Preliminary 
Measurements" and appendix A. 


The temperature differences obtained range from 8° F to about 110° F. 
Greater temperature differences are not desirable because the higher the 
value of At the greater the error introduced by the evaluation of the 
film temperature and the evaluation of the radiation, free convection, 
and conduction losses. The range of temperature difference was also 
limited by the method of calibration used for the heating unit. 


PRELIMINARY MEASUREMENTS 


Thermocouple Calibration 


The thermocouples used were calibrated against National Bureau of 
Standards standarized thermometers. When compared with the calibration 
tables from NBS (reference 28), the calibration for the copper-constantan 
thermocouple showed a discrepancy of 1 percent which is allowable 
according to standard practice. The copper-constantan thermocouple was 
made of one copper wire with two constantan wires in parallel so as to 
eliminate the effect of nonhomogeneity of the constantan wire (refer- 
ence 29). The thermocouple made from Cupron and Cromin D and used for 
particle-temperature measurements was made of only two wires since it 
was essential to maintain minimm wire cross section. 


Calibration of Heating Unit 


The calibration of the heating unit was made in the same manner as 
the test runs previously described except that the particle was enclosed 
in an evacuated glass bulb, the air in the column was stationary, and 
the cold junction of the thermocouple was placed in a melting ice bath 
in a thermobottle instead of suspended in the air stream. The glass bulb 
was suspended inside the glass column by the thermocouple wires. If the 
bulb were completely evacuated, the heat loss from the particle by 
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convection would be eliminated and the total loss would be due to radia- 
tion if conduction along the wires is negligible. When the induction 
coil was energized the particle temperature increased making it possible 
to obtain a plot of temperature rise against time. In addition to this 
heating curve, a cooling curve was obtained by observing the decrease in 
temperature with time when no energy was supplied to the heating coil. 
By combining these two curves the total rate of heat generation in the 
particle can be determined when the mass and the specific heat of the 
particle are known. The calibration of the induction coil was repeated 
for each size of particle used sinee the efficiency of the unit changes 
with the size of the particle. The calibration curves are shown in 
figure }. 


It is impracticable to obtain complete evacuation and maintain it 
for a period of time long enough to take a series of readings for this 
calibration. Therefore, in the actual calibration:a pressure of 7 milli- 
meters of mercury.was maintained. The heat loss due to free convectionl 
inside the glass bulb was then calculated from hr. at 1 atmosphere 


using the ratio of hr. at 0.01 atmosphere to hp, at 1.0 atmosphere 


equal to approximately 1/2. This value was based on the recent work of 
Gordon in which he studied the heat transfer by free convection from a 
piece of fine wire to the surrounding air at low pressures (reference 22). 
The value of hp. at 1 atmosphere was approximated in the present experi- 
mental apparatus with no air blown through the column and temperature 
differences not large enough to cause appreciable radiation. The heat 
loss by radiation was then obtained by deducting the heat loss due to 
free convection from the total loss which was obtained from the cooling 
curve. The heat loss due to radiation determined from these measurements 
agrees fairly well with the calculated value if a value of 0.7 is used 
for the emissivity of the steel particle in the usual Stefan-Boltzman 
equation : 


τι ME /τοὺ; 
αι = (0.7)(0.173) 755) - log (1) 


The above equation was then used for calculating the radiation corrections 
for the calibration of the heating unit. The conduction loss in the 
thermocouple wires was considered to be negligible (see appendix A). 

Even if this loss were. appreciable it would appear in both the calibra- 
tion runs and the test runs tending to eliminate this error. 


ithe free-convection effect is significant when temperature dif- 
ference is relatively large as in this case. 


aa amm ama “ὁ, À 


10 NACA TN 2867 


Estimation of Specific Heat 


From available data it can be seen that variation of the minor com- : 
ponents in carbon steel has negligible effect on the specific heat. A 
value of 0.115 Btu per pound per °F, based on the composition specified 
by one ball-manufacturing company, was obtained from the "Metals Handbook" 
(reference 30). However, since the compositions of all three particles 
were not available, it was desirable to have & simple calorimetric meas- 
urement of their specific heats in order to detect any significant dif- 
ference of composition. A simple water calorimeter made from a thin- 
wall test tube with an evacuated jacket was used. Twenty milliliters 
of distilled water were heated to approximately 2009 F and poured into 
the calorimeter. The steel particle was then submerged and moved up 
and down slowly by a mechanical device to provide agitation. A cooling 
curve for the water and a heating curve for the particle were obtained 
simultaneously by means of thermocouples. The specific heats determined 
in this manner for the various particles were in agreement within the 
experimental error and were close to 0.115. 


Pitot-Tube Measurement 


The velócity of the air stream at the center of the column was ^ 

measured by means of a pitot-tube probe using a micromanometer to detect 
. the pressure difference between the total- and static-pressure tubes. 

The velocity at the center of the column was also calculated from the 
average velocity obtained from & rotsmeter by using values of the ratio 
Umax/Uay at various flow conditions based on the data of Nikuradse 
(references 21 and 31) for the radial velocity distribution near the 
entrance of a column corrected for the presence of the particle at the 
center of the column (see reference 24). 


The calculated and experimental values, compared in figure 5, show 
good agreement, which justifies the use of calculated values especially 
at low velocities where the micromanometer fails to give accurate readings 
because of the extremely small pressure difference of the pitot tube. 

(See section entitled "Precision of Measurements." ) 


PRECISION OF MEASUREMENTS 


The accuracy of the data obtained from the measurements is affected 
by two types of error, systematic and random. Systematic errors were 
introduced in the instrumentation and its calibration, and random errors 
are due to inconsistent behavior of equipment and to human fallibility. 
While the systematic errors can be analyzed on a rational basis, the 
random errors can be detected only by the scattering obtained in the 2 
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correlation of data. Although every effort was made to reduce random 
errors to a minimum, a large number of runs were used in order to obtain 
compensation for this type of error. 


Temperature Measurement 


Since the cold junction of the thermocouple was placed in the air 
stream, the temperature difference between the particle and air stream 
was measured directly. The experimental accuracy is determined by the 
smallest graduation of the potentiometer reading. This gives a maximum 
error of 2.5 percent in the temperature readings. If an error of 1 per- 
cent is allowed for the calibration curve of the thermocouple used (see 
section entitled "Preliminary Measurements"), the total error will then 
be 3.5 percent for the temperature difference between the particle and 
the air stream. The effect of the location of the hot and cold junctions 
on the accuracy of the measurement is negligible because the temperature 
of the particle can be considered uniform throughout as shown in 
appendix B and the temperature gradient of the main stream was small 
compared with the temperature difference measured. 


In the measurement of the air-stream temperature at the outlet of 
the rotameter and at the test section the temperature readings were 
large enough to make the experimental error insignificant. 


Calibration of Heating Unit 


In the calibration of the heating unit the power input to the unit 
was determined from the plate voltage and current of the oscillator. 
The maximum error, when the power input was smallest, was 6 percent. 
This results in the same percentage error in the heat generated in the 
particle since the calibration curves relating the energy input to the 
heating unit and the heat generated in the particle show essentially 
constant slopes (fig. 4). However, the fact that the maximum error was 
not present in most cases can be shown by reproducibility tests. 


In computing the rate of heat generation Q in the particle, the 
slopes of the heating and cooling curve were used. Evaluation of these 
Slopes introduced a random error which cannot be readily analyzed. This 
error combined with the error due to instrument readings can be deter- 
mined from the maximum deviation of the data from the best curve drawn 
for the calibration. It is observed from figure 4 that the largest 
deviation of about 6 percent occurred in heating the 1/8-inch-diameter 
particle in the 3-inch-diameter column while in most other cases 2- or 
3-percent deviations are observed. The value of Q obtained from the 
calibration curve will then have a maximum error less than 6 percent. 
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Another random error which may enter into the value of the heat 
generated in the particle is that caused by inconsistency in the loca- 
tion of the particle during the test run and during the calibration. 
This inconsistency, if large, would change the efficiency of the heating 
unit and thus invalidate the calibration. Such errors were reduced to 
a minimum by keeping the magnetic field strength inside the induction 
coil as uniform as possible for a length of column several times the 
diameter of the largest particle. The effect of particle position was 
then tested and found to be less than 3-percent deviation for a dis- 
placement of 1/4 inch along the axis of the column and about 2-percent 
deviation for a horizontal displacement of 1/8 inch from the center. 
During the test runs the variation of position was much less than these 
values and, therefore, the error was well within the range i3 percent. 


The maximum error in the value of Q,.as well as in the measurement 
of temperature difference At, has a direct effect on the computation of 
the heat-transfer coefficient h and hence on the value of St. The 
combined error of these two quantities, Q and At, results in a maximum 


value of [1 + 0.05)/(1 - 0.035)| - 1 or 9 percent when the two errors 
are in opposite directions. 


Velocity Measurement 


The average velocity of the air stream and the maximum velocity at 
the center of the column were determined by different means. The maximum 
error in the average-velocity measurement determined from the least count 
of the rotameter scale is about 3 percent. The maximum velocity at the 
center of the column was measured by means of a pitot-tube sensing 
element connected to a micromanometer with a design sensitivity of 
0.00027 inch of water. Because of the relationship between the velocity 
and the pressure difference obtained from the manometer, this instrument 
will not give accurate measurement at velocities less than 2 feet per 
second. However, for air velocities greater than 2 feet per second the 
error becomes negligible except for error introduced by the inconsistent 
behavior of the apparatus. As shown in figure 5, this error is about 
3 percent. This error does not enter into the final result if Umax 18 


calculated from  Uay. 


The error in the velocity measurement affects both the value of St 
and Re but in a different manner. This effect can be determined by 
examining the relationship of these two dimensionless groups (fig. 6(a)) 
which can be approximated by the simple exponential equation 


St = Re (2). 


NACA TN 2867 | 13 


The effect of an error of 3 percent in the velocity measurement on the 
relationship given by equation (2) will be about 1 percent, since the 
value of M is approximately 0.6. When this error is added to the 
other possible errors discussed previously, a maximum total error of 
about 10 percent is obtained. 


ANALYSIS OF DATA 


Typical data obtained from this investigation are presented in 
table I. The heat-transfer coefficient was calculated from the general 
equation 


Q = hA(At) (3) 


Values of Q, the rate of heat transfer from the particle to the air 
stream for steady-state heat flow, were read from the calibration curves 
(fig. 4) of the heating unit and corrected for radiation loss by equa- 
tion (1). Values of At were measured directly by a thermocouple in 
the particle, with its cold junction in the surrounding air stream. The 
total surface area of the spheric particle is A. It is obvious that the 
heat-transfer coefficient obtained in this manner is an average coef- 
ficient over the entire surface of the particle. The heat-transfer 
factor (Stanton number St), was computed from the heat-transfer coef- 
ficient by means of the equation 


_ h 
ne (1) 


The heat-transfer factor was used in correlating the data with the values 
of the Reynolds number Re and also in the comparison of heat transfer 
with the skin-friction factor for the momentum-transfer phenomenon. 

The use of this dimensionless factor as a means of correlating heat- 
transfer data was suggested by McAdams and Drexel (reference 32), 
although most available data were presented in terms of the Nusselt 
number Nu. 


The properties of the air were evaluated at the film temperature 


which was assumed to be the arithmetic average of the particle tempera- 
ture and the temperature of the air stream at the test section, or 


Tp = Tair + Z(At) (5) 
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The values of the properties of the air were based on those given in 
the NBS-NACA "Tables of Thermal Properties of Gases" (references 33 
and 34). A correction was applied for the moisture content in the air 
determined by taking a humidity reading at various times during the 
series of test runs. 


SKIN FRICTION ON SURFACE OF A SPHERE 


Since momentum transfer along a surface is the result of skin 
friction, an evaluation of this skin friction can be used in attempting 
a correlation between momentum transfer and heat transfer. The direct 
measurement of skin friction on the surface of & sphere with & stream 
flow in the Reynolds number Re range of 50 to 1000 was impracticable 
because of the minute apparatus required for ihe small-size sphere used. 
Furthermore, since the validity of experimental data for skin friction 
has been questioned (reference 17) it was decided to calculate values of 
the skin friction for use in the proposed correlation. 


For the boundary layer the Navier-Stokes equation can be simplified 
(reference 21) into the form 


2 
αϑα δα iae, V A (6) 


Ox ày "px 3 


Where u and v are velocities along two perpendicular axes x and y 
and the plane of (x,y) is the plane of motion; p and "V are the 
density and kinetic viscosity of the fluid. Since the total change of 
pressure along the normal to the surface throughout the boundary layer 


is negligible,? 


— 
— 


ΟΡ aU 
a (7) 


_i 
p 


where U is the velocity in the main stream just outside the boundary 
layer. 

"This is also generally true for a curved surface. However, when 
the boundary-layer thickness is fairly great, the pressure gradient 
across the layer that is required to balance the centrifugal force 
produces an appreciable pressure drop between the outside of the layer 
and the surface. The fundamental assumption of the boundary-layer 
theory will be, to this extent, violated. 
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These equations combined with the continuity equation 


du , 9v _ 
ὃς MS = 0 (8) 


provide the basis for calculation for all characteristics of the boundary 
layer. While equations (6) and (7) were developed for two-dimensional 
flow, they can also be used for flow around the surface of a solid of 
revolution such as a sphere (see references 21 and 35). In this case, 
the continuity equation should have the form 


9(rou) , (rov) BR, 


I = (9) 


where το is the perpendicular distance from the axis of revolution to 


the surface. As v is small compared with u, equation (9) can be 
expressed approximately as 


à 
Sin Z-0 (10) 
o 


Approximate Solutions 
Millikan has applied the momentum equation, which can be obtained 
by integrating equation (6) with respect to y between O and δ, to 


spheres and has derived, after certain simplifications, the equation 
(reference 35) 


ὃ I 2 9 1 dro [ 
pu^ dy - U — pu dy - — ——|U pu dy - 
ox Jo 21 To 9x V Jo 


2 2 op 2») 
put dy} = -8 .pn[ sz (11) 
[ ο òy) y-o 


where x is measured from the forward stagnation point along a meridian 
curve as shown in figure 7. 


ge gen, 
MM —— u et ep m m e. P a e ne 


16 NACA TN 2867 


This equation for a sphere differs from the momentum equation for 


dr 
two-dimensional flow only by the terms containing = =, Hence the 
ο 


same equations used in the Pohlhausen-Kärmän method (references 36 

and 37) can be applied to the sphere if these terms are added. In the 
present calculation of skin friction a further modification by Holstein 
and Bohlen (reference 38) was used. 


Velocity Distribution 


In the calculation of the skin friction the velocity distribution 


outside the boundary layer is required as shown by the term + 9. which 


is equal to us. This velocity distribution must be experimentally 


determined before the calculation can be performed. Such data for spheres 
with values of Re in the range of 50 to 1000 are not available.  How- 
ever, as mentioned by Tomotika (reference 39), as long as the velocity 

of the undisturbed stream is below the critical value for the sphere used, 
the velocity distribution derived from one particular pressure-distribution 
curve experimentally determined for one stream velocity may be taken as 
representative of the velocity distribution around the sphere. Asa 
verification of this statement, the data available in the literature for 
the pressure distribution surrounding a sphere and a circular cylinder 
(references 19, 21, and 40) are shown in figure 8. The pressure- 
distribution ‘curves for cylinders are available for values of Re as 

low as 2800. By observation of these curves two general trends may be 
noted. First, although the pressure-distribution curves for various 
Reynolds numbers show deviations at the surface near the maximum velocity 
point (a minimum pressure shown in the figure), at the downstream surface 
all the curves tend to approach a constant pressure which is not greatly 
affected by the value of Re. Second, the smaller the value of Re, 

the sooner the pressure curve recovers from the minimum values resulting 
in a longer flat portion of the curve. The curves for low values of Re 
have only one minimum point if the minor irregularities, such as seen in 
curve (II), are neglected. Considering these trends, the shape of the 
curve for even smaller values of Re would be expected to be similar to 
the shape of curve (I). Furthermore, calculations of skin friction based 
on curves (I) and (II) show a deviation of only 3 percent for the total 
skin friction on the surface of a sphere. It is, therefore, considered 
justifiable to use curve (I) as the pressure-distribution curve for 
calculation of skin friction in the Re range of 50 to 1000. 
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Result of Calculation 


The skin friction on the surface of a sphere can then be calculated 
for each value of x from the stagnation point to the point of separa- 
tion where the boundary layer leaves the surface. The point of separa- 
tion, 0 = 86°, was determined from the calculation and is also indicated 
by the pressure-distribution curve. Beyond this point, the local skin 
friction is small as shown by Geidt (reference 19) and in reference 21. 
Because of the existence of vortex rings which become unstable (refer- 
ences ll and 12) at values of Re above 150 and the fact that well-defined 
turbulence does not exist until much higher values of Re are attained, 
it is probable that the usual concept of stagnant gas films for turbulent 
flow is not strictly applicable to this range of Re. As an estimation, 
2 percent was added to the skin friction calculated for the laminar 
boundary layer. This value is the same as that determined for circular 
cylinders by Thom (reference 21). 


For the correlation of skin friction with the heat-transfer data, 
a dimensionless factor will be used. This factor, designated as Cf, 
is defined as 


T dÀ 
Total frictional force on surface H 
Ce/2 = ee Dee on eee ὦ (12) 


ρῦ ο” Α Όσο Αρ 


where A is the area, τ is the shearing stress on the surface, or 


(a) ; and Ug is the velocity of the undisturbed stream. The area 
y=0 


in this equation is the total surface area of the sphere and is, there- 
fore, different from the projected area on which the drag coefficient 
ls generally based. Furthermore, the total frictional force has to be 
distinguished from the friction drag. The former is the summation of 
the local frictional forces (shearing force) over the whole surface, 
while the latter is the resultant of the frictional forces in the direc- 
tion of flow of the main stream. In order to compare with the heat- 
transfer process on the same surface, the skin-friction factor Cr is, 
therefore, based on the total frictional force. The direct correlation 
between the heat-transfer factor and the friction drag, as suggested by 
previous investigators (references 20 and 41), is questionable. 


The results of the point-by-point calculation of the values of 


| L5 Js ‘were integrated graphically over the area covered by the 
ρῦο 


ee eee 
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calculation, and then divided by the total surface area of the sphere, 
giving the value l 


P E ls ee : fis n 
——zj Re dA 21 e|r^ sin 0 dé 
O E | _ OP EE ) - 0.63 


A are 


Adding 5 percent for the portion of the surface beyond the separation 
point, this becomes 


A 
RESE 
O Ex , (=!) (e)? =- 0.66 


A 


which, after rearranging, results in the final expression 
Ce/2 = 0.66 (Re) (13) 


DISCUSSION 
Comparison of Experimental Results with Previous 


Heat-Transfer Data 


The results obtained from the present investigation are presented 
in figure 6(a) as a plot of St against Re and in figure 6(b) as a 
plot of Nu against Re. The curve in figure 6(a) is represented by 
the empirical equation 


St - (3.10/Re) 4 |ο.ογ/(πογθ-3 (14) 


with a maximum deviation from the curve of 5 percent which is within the 
experimental error. For comparison the results from previous investiga- 
tions on spheres are also shown in figures 6(a) and 6(b). 


Williams has suggested & curve for spheres based on several sets of 
data reported in the literature (reference 3).2 As previously mentioned, 
the ‘application of this curve in the range of Re from 50 to 1000 is 
limited. 


JThe curve given in reference 3 was corrected according to a private 
communication (1950) with the author. 


————€€——————— Se + os —— -- —-- --- 
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Kramers has made a study similar to the present one, using air as 
one of the mediums flowing around a sphere (reference 24). His investiga- 
tion covered & Re range from 200 to 2000 and gave a curve having a 
slope parallel to the present results as seen in figures 6(a) and 6(b). 
He found a linear relationship between Nu and (Re)Ü-? which may be 
represented by the equation 


Nu = 3.2 + 0.59(Re)9-? (15) 


However, this equation gives consistently higher values than the present 
results with a deviation of 50 to 55 percent. 


For comparison, equation (14), which represents the present data, 
can be converted to a form similar to equation (15) by using a constant 
value of 0.67 for Pr corrected for moisture in the ambient air. This 
gives the equation 


Nu = 2.1 + 0.37(Re)9-? l (16) 


It is evident that at low values of Re equation (16) gives a smaller 
value for Nu than does equation (15), and the value is closer to 2 
which has been shown (references 1 and 23) to be the limiting value 

of Nu as Re becomes very small. Since deviations from the theoret- 
ical value are generally considered to be due to free convection this 
also indicates that the effect of free convection was within the experi- 
mental error. 


Ranz and Marshall (reference 42) have recently suggested the 
equation i 


0.33 


Nu = 2.0 + 0.60(Re)" > (Pr) (17) 


to correlate the heat-transfer data for the system of a water drop in 
air. For dry air where Pr = 0.71, equation (17) becomes 


Nu = 2.0 + 0.53(Re)°°? (18) 


This can also be converted to a form comparable with equation (14) to 
give 


St = (2.82/Re) + πο. (19) 


en 
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These equations give St and Nu curves, as plotted in figures 6(a) 
and 6(b), parallel to the curve obtained from the present investigation, 
but show a maximum deviation of about 28 percent in the range of Re 
studied. No satisfactory explanation is offered for this deviation. 
However, as pointed out by Ranz and Marshall (reference 42), the relia- 
bility of this correlation depends to a great extent on the values used 
for the transport properties of the fluid. The values of thermal con- 
ductivity of dry air used in the present calculations were obtained from 
the NBS-NACA "Tables of Thermal Properties of Gases" (references 33 

and 34). These values are about 3 percent higher than those used by 
Ranz and Marshall (reference 42). 


The theoretical curve shown in figures 6(a), 6(b), and 9 was 
obtained by Drake, Sauer, and Schaaf (reference 2) by solving the steady- 
state energy-flow equation for spheres.  Expressed in spherical coordi- 
nates, as illustrated in figure 7, this equation is 


ea 


oT ı OF Y Or _ 
"RS, * M9 Fo "rsinjod 


1 3[25|, faino E|- 1. A 
JE ἫΝ E) sin 0 Sl =) r? sin? 0 e on 


In obtaining this solution, several simplifying assumptions were made 
among which is the assumption that the velocity of the fluid around the 
sphere is everywhere parallel to the surface and equal to the undisturbed 
velocity in the main stream, or Up = Ug. 


It is evident that this approximate solution deviates from the 
true solution of the energy equation when the value of Re is large 
enough that this assumption is no longer valid because of the fact that 
the fluid stream separates from the particle surface. However, at small 
values of Re where the velocity term does not have an appreciable effect 
on the result, it would be expected that this approximation will approach 
the true solution and, therefore, the curve will represent the true 
values at low values of Re. With this in mind, it is expected that the 
experimental curve will approach Drake's curve as Re becomes small. 
This is seen to be the case. This agreement also leads to the possibility 
of extending the present data below Re = 50. 


Another analytical solution, based on the thermal-boundary-layer 
theory, by Kudryashev (reference 23) is 
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Nu = 2.0 + 0.33(Re)9°? (21) 


This equation gives still lower values for Nu than those obtained from 
either equation (16) or (18). 


Comparison of Skin Friction with Heat-Transfer Data 


The existence of an analogy between heat transfer and momentum 
transfer on the interfacial surface of a solid and a flowing fluid makes 
it possible to correlate these two transfer phenomena. The relationships 
that’ have been found useful in the cases of fluid flow in conduits and 
on flat surfaces suggest the desirability of such a correlation for the 
case of spherical particles. In the case of fluid flowing around a 
blunt object, the presence of impact pressure on the surface of the 
object also contributes to the momentum transfer, and therefore must be 
distinguished from that momentum transfer due to skin friction. The 
term "momentum transfer" as used here refers to that part of the total 
transfer which is due to skin friction alone since it is this momentum 
transfer which is analogous to heat transfer as shown in appendix C. 


The factors which are comparable in these two transfer phenomena 
are Cr /2 and St. Since the flowing medium used in this investigation 
is air, the Prandtl number Pr is close to 1 and the ratio St/(Cr/2) 


can be shown to be unity. 


As shown in figure 9, the heat-transfer factor St curve resulting 
from the data obtained in this investigation lies between the theoretical 
curve calculated by Drake, Sauer, and Schaaf (reference 2) and the Cr/2 
curve calculated in the section entitled "Skin Friction on Surface of a 
Sphere." Over most of the Re range studied the experimental curve 
shows better agreement with the friction-factor curve than it does with 
Drake's theoretical curve. The deviation between the Cf/2 curve and 
the St curve at low values of Re is explained by the increase of the 
boundary-layer thickness. The assumption of negligible pressure drop 
across the layer, used in the calculation of the skin friction, becomes 
invalid as the boundary layer becomes thicker. The fact that the friction 

IN 


ccording to Von Kärmän's equation (references 11 and 12), 


Cr/2 


14-5 Gre tex -ᾱ1 loge [i + (5/6) (£r. - 1] 


then for values of Pr between 0.71 and 1 


St = 


St = Cp/2 


m nn ——— — 
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factor at low values of Re should be higher than the values calculated 
by the boundary-layer theory is indicated by the curve obtained from 
Stokes' equation for very low values of Re. 


it is of interest to note the temperature profile around an 
evaporating water drop as reported in the recent work of Ranz and 
Marshall (reference 42) who point out that the thickness of the boundary 
layer varies with the position on the surface of the particle and is 
comparable with the diameter of the drop. This would seem to invalidate 
the calculation of skin friction based on the assumption of a thin 
boundary layer. However, this variation of boundary-layer thickness 
of the streamlines from the particle surface and, therefore, the calcu- 
lation of skin friction from the stagnation point to the point of stagna- 
tion as used in the section entitled "Skin Friction on Surface of a 
Sphere" remains satisfactory. 


It may be pointed out that the absence of an accurate method for 
determining the fluid dynamics in the range of Re studied has prevented 
obtaining better agreement in the comparison of the momentum-transfer 
with the heat-transfer factor. 


The above considerations seem to indicate the existence of the 
relationship St = Cr/2 as expected from theoretical analysis (see 
appendix C). This is also in agreement with Sherwood's findings in 
the case of a cylinder (reference 20). In both cases the relationship 
applies for gases where Pr is close to unity. The existence of this 
simple relationship between heat and momentum transfer for a spherical 
particle in a gas stream discloses an additional application of the 
theory of analogy between these two transfer phenomena. It should be 
expected that this analogy would apply to any blunt object having a 
smooth surface in addition to its application to streamline bodies and 
flat or curved surfaces. 


Although there is no apparent analogy between the drag coefficient 
Cp and the heat-transfer factor St inthe Re range studied, they 
both show similar behavior in varying with Re. This is believed to 
be due to the fact that the form drag does not have a significant effect 
on the variation of total drag with change in Re since the air-stream 
velocity in this Re range is far below the critical value and, there- 
fore, turbulence at the wake of the sphere is not an essential factor. 
The total-drag coefficient thus varies with Re in approximately the 
same manner as the surface friction drag. Since the total-drag coef- 
ficient is measurable and is available in the literature (references 21, 
43, and 44), an empirical relationship between Cp and St, based on 
the present data, was developed (see fig. 10) resulting in the equation 


St = (Cp/8) - fo.10/ (πω) 13] (22) 


¢ 
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The maximum deviation of this equation from the St curve obtained from 
the present experimental data is about 8 percent which is within the 
experimental error. Such an empirical correlation between the heat- 
transfer factor and total-drag coefficient may be useful for the deter- 
mination of heat-transfer rates from drag measurements especially in the 
case of submerged bodies which are not true spheres. Such application 
may be valid where a correction factor (or form factor) is employed in 
the above equation to allow for nonsphericity. 


CONCLUSIONS 


The conclusions which can be drawn concerning heat and momentum 
transfer for the system of a spherical particle in a moving air stream, 
based on the results obtained in this investigation, are: 


l. A correlation of the heat-transfer factor St against the 
Reynolds number based on particle diameter Re was obtained for the 
range of Re from 50 to 1000. This correlation is represented by the 
equation 


St = (3.10/Re) + [ο.55/ (3909-21 


At the extremes of the limits for the Re range studied the experi- 
mental data are in agreement with values predicted by theory which 
tends to substantiate the validity of the experimental curve. 


2. Extension of the experimental curve into the region below values 
of Re of-50 by means of Drake's theoretical equation is indicated. 


3. The satisfactory agreement between the calculated values of the 
momentum-transfer factor Cg/2, where Cp is the skin-friction factor, 
and the experimental data at the higher values of Re studied suggests 
the relationship St = Cg/2 for the case of air flowing around a sphere. 


4. For practical application in the Re range of 50 to 1000, the 
empirical equation 


St = (Cp/8) - ο 15] 


relating the heat-transfer factor to the readily measurable total-drag 


coefficient Cp may be useful for particles deviating slightly from 
true spherical shape. 


University of Florida 
Gainesville, Fla., December 10, 1951 
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APPENDIX A 
HEAT LOSS THROUGH SUSPENSION WIRES 


The conduction losses along the thermocouple wires were minimized 
by using extremely fine wires. However, for precision the magnitude of 
these losses is considered. The losses due to conduction to the ends 
of the wires include losses from the shellac-coated wire surface by 
convection and radiation. Thus the maximum loss by conduction would be 


Qc = -kAg(aT/aL); o (A1) 
The Fourier equation for conduction (reference 15) can be written as 
OE are 


where the term b°T represents the rate of heat loss through the surface 


element to the surrounding medium, b^ is a constant to be determined, 
€ is time, and T is the temperature of the wires referred to that of 
the surrounding medium. For steady-state heat flow, this reduces to 


(A2) 


c DELE | (A3) 


. Equation (A3) can be solved by substituting ο for T. This 
results in the solution 


v = Be>L/VE + qe-bL/Va (Ah) 


where B and C are constants. The boundary conditions for the longi- 
tudinal temperature distribution along the wires are Το at L= 
and T=At at L=0O. 


Therefore, B=0, C = At, and 


T = (at) (eL Va) (45) 


from which 


(ar/àL), ο = -(At)b/ve (AG) 
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In order to find the value of b, consider a finite section of length AL 


of the wires. A heat balance can be written as 


kAg AL = hHTgPAL . (AT) 
dL 
or 
d?T 1Ρ 


where h is the sum of radiation and convected film conductance which 
is assumed to be constant over the range of temperature considered, P 
and Ag are the perimeter and the cross-section area, respectively, 
and Tg is the temperature at the surface of the coated wires. By 


combining equations (A3) and (A8), b^ can be expressed as 


hP 
be = in. (Is/T) (49) 


The value of Tg/T determined from a heat balance across the shellac 
film is about 0.5, assuming a value of 15 for h at the surface of the 
shellac film (based on the data of Gordon (reference 22)). By sub- 
stituting the proper values of P, Ag, h, and k, equation (A9) gives 


b/Vva = 90 ft 


Then 


and 
Qe = kAg(90) (^t) 


Under the test condition such that the minimum amount of heat 

(0.030 Btu/hr) is transferred corresponding to a value of At of 10° F, 
the conduction loss will have the maximum effect. This maximum error 

in the value of Q due to conduction loss may reach 3 percent. In most 
cases the amount of heat lost by conduction was found to be negligible. 


-- — πι -- νι — M —— ος — — — y aa ee — ————————MS ti RR mc g κ 
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APPENDIX B 


TEMPERATURE INSIDE STEEL PARTICLE 


Temperature Gradient during Heating 


It is well-known that when a metallic particle is heated by an 
induction coil, the eddy current induced in the particle is concentrated 
on the surface and the interior of the metal receives heat from the skin 
by thermal conduction. Because of the high conductivity of steel the 
heat per unit volume generated in any concentric spherical element is 
assumed to be the same as for any other element. Thus, 


Q/R? = Qo/ Ro” (81) 


where Q, is the total heat generated per unit time in the sphere with 
a radius of R, and Q is the heat generated per unit time in & con- 


centric sphere with radius of R, which is smaller than Ro. If Αρ 
is the surface area of the sphere and A is that of any concentric 
sphere with a radius of R, then 


B= ale) RE 68) = 


From Fourier's law for steady conduction 


dt _oRı (B3) 
dR — kgA Ay Ro kg 
Then 
R 
Qof 1 
dt = - ---ί------ R dR 
Ε κοίτης) i 
or 
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where t, is the surface temperature of the spherical particle and Κα 
is the thermal conductivity of the particle. At equilibrium, when a 
steady state is established, the heat generated is equal to that trans- 
ferred to the ambient air stream : 


Qo/Ag = h(At) 


where At is the temperature difference between the particle and the 
main air stream. Substituting this value in equation (BÀ) gives 


2 2 
E n{Ro~ - R : 
Ὁ - tg = oo 2 RE) (85) 
Let t. be the temperature at the center of the sphere. Then 
hR hR„.\/k 
ο O f 
tos tys (δύ)σε- = (ze) (B6) 


where kp is the thermal conductivity of the ambient fluid. It is 
evident that kr for air is much smaller than kg (kr/kg is approxi- 
mately 1/2000) and since the value of Nu which is equal to 2hRo/kp 
is never more than 20 for the range of Re studied, the value of 

te - tg is less than 0.25 percent of At.? Hence the error introduced 
by considering tc equal to tg is negligible for the values of At 
encountered in this investigation. 


Temperature Gradient during Cooling 


In calibrating the heating unit, cooling curves of the particle as 


well as heating curves were used. The fact that the temperature measured 


by attaching a thermocouple at the center of the particle can be used in 
place of the surface temperature must be established. 


To determine the possible error introduced by the assumption of a 
uniform temperature within the particle at any instant the temperature 
of the rarefied gas surrounding the particle inside the evacuated bulb 
is assumed to be uniform and equal to the surface temperature tg. If 


2Since the temperature at the surface equals the temperature at the 
center & uniform temperature must exist throughout the particle which 
justifies the original assumption made for the heat distribution. 


mem —— ——— — —M— mn ————— ———— nn — Ha ——— MM — M M P ——— MM € —MM Màn er P σσ li Ange PUPA —— —grÀ M MÀ 
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ts is suddenly reduced to some temperature lower than the initial 
temperature of the particle” there will be an instantaneous temperature 
gradient established within the particle because of the thermal resist- 
ance to the conduction of heat from the center to the surface of the 
particle. Then, from the Fourier equation of conduction (reference 45), 


t+ - ts 
and 
ta -t - T 
c = Seb + Pele + pe. l .) = Bale) (B8) 


where t. is the temperature at the center of the particle, t4 is the 
initial uniform temperature of the particle, tg is the average tempera- 
ture of the particle at time e, and g is defined by 


g- 1^a.c/R,* 
If the time interval e is 2 seconds and the thermal diffusivity a 


of steel is 0.570 square foot per hour, then for the largest spherical 
particle used (Βο = 5/16 in.) 


g = 4.61 


From the table of values for the functions B(g) and Balg), (refer- 
ence 45), it is found that 


B(g) = 0.020 
Ba(g) = 0.006 


Then, by combining equations (B7) and (B8), 


(ta - tg) - (tc - tg) = (0.020 - 0.006)(t4 - tg) 


Thus the error in this case which would be the maximum error encountered 
in the tests will be 1.4 percent of the temperature change from the 
initial temperature to the surface temperature. 


his is a hypothetical condition which represents the most unfavor- 


able condition resulting in maximum error. 
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APPENDIX C 
HEAT- AND MOMENTUM-TRANSFER ANALOGY 


Considering a laminar boundary layer on the surface of a sphere 
where no fluctuations in velocity and temperature exist, the analogy 
between heat and momentum transfer (references ll and 12) can be seen 
from the similarity in form of the following equations. For heat transfer 


dt 
α/ορρ = -Q dy (C1) 
and for momentum transfer 
r/o =v ἂν (c2) 


where y is in a direction perpendicular to the surface. 


These two equations will be directly proportional if V and a are 
numerically equal, or 


v/a = ucp/k = Pr = 1 


which is approximately satisfied by a gas. If in this layer q and τ 
are assumed to vary with y ina similar manner, the velocity and tem- 
perature at the boundary are Ug and tm, respectively, and t' is the 
temperature of the fluid referred to a surface temperature of zero, then 


tm : 
a | dt = | dy (C3) 
o (α/ορρ) ο 


ο , 
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Therefore, 
ὕο/(τ/ϱ) = epe(At)/a (C5) 


By using the definitions of the skin-friction factor Cpr =T / (puc*/ 2) 
and heat-transfer factor St = Ἠ/ορρῦο = q/cppUQ(At), the above relation- 
ship can be transformed into the familiar form 


Cp/2 = St (C6) 


This relationship was derived on the assumption that the boundary 
layer for the velocity field is the same as that for the thermal field. 
This is true, however, only for gases. In general the hydrodynamic 
boundary layer and the thermal boundary layer have different thicknesses, 
5 and δὲ, respectively. For the determination of the thermal-boundary- 
layer thickness and the heat transfer through this layer, a heat-flow 
equation for the thermal boundary layer can be used. The desired equa- 
tion is derived from a heat balance for a volume element of the fluid at 
the surface of the sphere, similar to the derivation of the momentum 
equation for the hydrodynamic boundary layer: 


1 l 
PC &l tfu dy - pc, &l t'u dy - xS) o -0 (CT) 


where tp is the temperature of the undisturbed stream and 1, the 


length of the volume element in the y-direction, is assumed to be greater 
than either ὃ or 64. By introducing the thermal diffusivity of the 
fluid a = k/cyp, equation (C7) becomes 


ὖ 
d dt! 
2 tpe t! = q| = c8 


By assuming the temperature profile in terms of y/5; to be analogous 
to the velocity profile in terms of y/5 and (8/5) « 1, the value 
of 94/5 can be obtained (reference +): 


δι/5δ = 1/ (2x)?/3 (C9) 


NACA TN 2867 31 


Since Pr for gases is smaller than 1, 6/5 becomes greater than 1, 

and the assumption made in obtaining equation (C9) is not strictly true. 
However, the error introduced by this approximation when the Prandtl 
number is close to 1 is very small, giving a value for δΙ/δ of approxi- 
mately l.l. It is, therefore, considered satisfactory to assume δε 
equal to 5 in the analogy equation, with a possible error of 10 percent. 


32 
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Rotamater reading, ou ft/min 
Thermocouple reading (I), αν 

Air temperature at rotameter outlet, OF 
Thermocouple reading (II), mv 

Air temperatura at test section, °F 


Air velocity at particle cross 
section, fi/gac 


Power input to heating unit, volt amperes 


Rate of heat transfer, Q, Btu/hr 


Thermocouple reading (III), mv 


Temperature differance, Ab, OF 
Film temperature, OF 
Reynolds number, Re 


Heat-transfer oosfficient h, 


Btu/tt* /nz /OF 
Hent-transfer factor, Bt 


Russelt number, Hu 


lonly typical rune are listed, The complete table of data is avnilable in the unpublished oopy of the annual report to KACA from the Univ. of Flarida 
submitted under contract MAw-6010 and entitled: "A Special Analytical Study of Heat Transfer Phenomena with Clouds of Small Particles in Air at Elevated 
B, Part I: Study of Heat Transfer with Single Particles" by J. M. Duncan, Y. 8. Tang, C. R. Wallace, and C. M. Oktay. This report is available 
for loan or reference in the Division of Hessarch Information, National Advisory Committee for Aeronautics, Washington, D. C. 


Designations: A - particle size, 1/8 in.; B - particle size, 1/8 in.; C - particle size, 5/8 in. Run mmber in 200 series representa P-in. I.D. of 
glass column used and run number in 300 series reprasents j-in. I.D. of glass colum used. SEA 
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TABLE I - Concluded 


1983 NL VOW 


section, ft/sec 


Power input to heating unit, 
volt ampares 


Rate of heat transfer, Q, 
Btu/hr 


i gi ον coefficient, h, 


Btu/ft"/hr 
Heat-transfer factor, 8t 
Husselt number, Nu 


Designations: A - particle size, 1/8 in.; B - particle size, l/h in.) O - particle size, 5/0 in. Run mwber in 200 series represents 2-in. I.D. 
of glass column used and run mmber in 300 series representa 3-in. I.D. of glass column used. SA” 
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Figure 1.- Apparatus for heat-transfer measurements. 
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Figure 2.- Schematic diagram of apparatus for heat-transfer measurements. 
1, blower; 2, surge tank; 3, air heater and temperature regulator; 
4, rotameter; 5, manometer; 6, glass column; 7, heating coil; 
8, Westinghouse R-F generator; and 9, potentiometer. 
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Figure 3.- Sectional view of spherical particle with hot junction of 
thermocouple attached. Four times full scale. 
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Figure 4.- Calibration curves for heating unit. 
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Figure 5.- Comparison of observed Umax and values of Umax 
from flow-rate measurement at center of column. 
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Figure 6.- Comparison of experimental results with available 
heat-transfer data. 
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Figure 6.- Conciuded. 
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(a) Orthogonal coordinates. 


(b) Spheric coordinates. 


Figure 7.- Orthogonal and spheric coordinates for spheres. 
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(VI) Reference 21 


(III) [Reference 21 


(VII) Reference 21 


Figure 8.- Pressure-distribution curves over a sphere and a cylinder. 
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Figure 9.- Plot of heat-transfer factor St or momentum-transfer 


factor Cpf2 against Reynolds number Re. 
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Figure 10.- Comparison of total-drag-coefficient function σηγθ with 
momentum-transfer factor C,/2 or heat-transfer factor St. 


NACA-Langley - 9-2-54 - 120 


